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Abstract 

We calculate production rates for massless (m, d) and massive (s, c, b) quarks in pure 
glue and quark gluon plasmas to leading order in the strong coupling constant g. The 
leading contribution comes from gluon decay into qq pairs, using a thermal gluon prop- 
agator with finite thermal mass and damping rate. The rate behaves as a|(ln l/a;5)^T^ 
when m,as —>■ and depends linearly on the transverse gluon damping rate for all values 
of the quark mass m. The light quark {u, d, s) chemical equilibration time is approx- 
imately 10-100 for g =2-3, so that quarks are likely to remain far from chemical 
equilibrium in ultrarelativistic nuclear collisions. 
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Ultrarelativistic nuclear collisions may present a unique opportunity to study the high 
temperature behaviour of hadronic matter in the laboratory. It is possible that in these 
experiments, the produced energy densities might be high enough and last for long 
enough for bulk, locally equilibrated, deconfined matter, or quark gluon plasma (QGP), 
to be created and studied. However, recent work has shown that, while the hot 
matter should thermalize in about 0.3 fm/c, it may be largely a gluon plasma (GP), 
with few quarks. Chemical equilibration is expected to take much longer than thermal 
equilibration, if it occurs at all. 

One of the proposed signals for QGP is copious production of strange particles 0, as 
the strangeness production rate is much higher in QGP than in normal hadronic matter. 
Similarly, it has also been argued that thermal charm production in a QGP could be 
large [Q. However, these rates have never been calculated using thermal field theory. 

In this paper, we calculate the thermal quark production rate to leading order in as 
for GP and for QGP. This gives the chemical equilibration rates for both the massless 
(u, d) and massive (s, c, h) quark components of QGP. We find that these rates are 
small; the light quarks (u, c?, s) chemical equilibration time is approximately 10-100 
for g =2-3, so that quarks are likely to remain far from chemical equilibrium in 
ultrarelativistic nuclear collisions. 

We take the idealized situation of a GP or QGP in thermal equilibrium. If such a 
plasma does exist, the massless gluons evolve into quasi-particles with effective masses 
of order gT [||. Being massive, these quasi-gluons decay into qq pairs. This situation is 
similar to the plasmon decay into vv pairs that is the dominant cooling mechanism for 
high density stars [If. 

For gauge field theories at high temperature, there exists a resummation program 
developed by Braaten and Pisarski that describes these quasi-particles. It is easy to 
check that this resummation program is almost the same for the cases of QGP and GP. 
The difference is that the quark contribution inside each hard thermal loop is omitted 
in the case of a GP, because there are no thermal quarks. For instance, the plasma 
frequency is 



for SU{N) gauge theory at temperature T, where g is the strong coupling constant and 
Nf is the number of massless fermion flavours. Thermal effects in QGP and GP are thus 
identical, except that Nf = for GP. 

The discussion of the different QCD corrections to the production of qq pairs is best 
illustrated by using "cut" diagrams The lowest order diagram is shown in Fig. 1. 

To disentangle the contributions from different graphs and from different discontinuities 
within the same graph, we consider only the leading term in the strong coupling constant 
g as g ^ 0. 

Following the usual treatment of Braaten and Pisarski ^ in such a situation, we 
consider the two cases when the gluon momentum is soft, Q ~ gT, and hard, Q ^ T. If 
it is soft, one has to use the effective quark-gluon vertices and gluon propagators which 
resum the hard thermal loops . One is therefore left with a quite involved calculation as 
several discontinuities can be taken in such a case. In principle, this situation should be 
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very similar to previous calculations such as photon emission from a QGP ||T0[- Indeed, 
the naive lowest order diagrams for the qq production are the two-gluon fusion processes 
shown in Fig. 2. However, it is well known that these diagrams have a logarithmic 
singularity when the quark mass, m, goes to zero [|n|]. The cross section for this process 
is 

^99^qq = ^^'^^ ^^^"^ S > m . (2) 

The divergence being of kinematic origin, it is clear that in the case of unthermalized 
and zero mass quarks, this singularity can be screened only by the thermal gluon mass, 
of order gT. As typical gluon energies are of order T, the quark production rate must 
then behave as g'^T'^hil/ g. When the unthermalized outgoing quarks are massive, the 
cross section for the processes in Fig. 2 is no longer singular and leads to a contribution 
of Oi^g'^T^) to the production rate. 

A consistent calculation would require a division of the processes into two pieces by 
the introduction of an intermediate cutoff gT <^ q* <^ T on the gluon line (as in the 
photon emission case). Then the hard part, calculated from the first two graphs shown 
in Fig. 2, gives a contribution of order gf^T^ ln(T/g*). On the other hand, the soft part, 
calculated from the discontinuity of the graph in Fig. 1 together with effective vertices 
and gluon propagator, must give g'^T'^ln{q* /ujq) with the same coefficient in front so 
that, by adding the two contributions, the q* dependence cancels. To be honest, we have 
not checked that the coefficient in front of the logarithm is the same for the soft and the 
hard pieces but we do not see any reason why it should differ. 

If the gluon momentum is hard and if one is using a bare gluon propagator, the rate 
vanishes for obvious kinematic reasons. However, as in the 7 z/z/ process, one could 
also use a resummed propagator for the hard gluon, leading to a finite answer (in fact 
a contribution of order g'^T^ as we shall see). But unlike the photon, the gluon has an 
anomalously large damping rate and if the hard thermal gluon mass starts to be relevant, 
one must also include the damping rate inside the hard gluon propagator. We show that 
the rate behaves then as g'^T'^ {In 1 / g)"^ for massless quarks and as g^T^ In 1/(7 for massive 
quarks. It is therefore larger than the soft gluon contribution and than the naive lowest 
order diagrams shown in Fig. 2. The precise reasons for this unusual behaviour will 
become clear as the calculation proceeds. 

The gluon propagator that we use is therefore 

g2-ReHT(L)(c^,g) + 22gV(L)' 

where P^^l) transverse (longitudinal) projector and a and b are colour indices. 

Here (and for the rest of this paper) the subscript T{L) refers to transversely (longitudi- 
nally) polarized gluons. We have included the thermal gluon mass at hard momentum, 
ReHT(L)(ci;, g), and also the anomalously large gluon damping rate, 7t(l) ]T^ - 

Using the cutting rules of Kobes and Semenoff , the quark (antiquark) production 
rate due to transverse gluon decay is given by 

R^=^ = g^ J^J^^ {2.f 5[K^ - m^] 6[{Q - Kf - m\ns{u) + 9{-q^)) 
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Here K and Q are the quark and gluon four-momenta respectively, m is the quark 
mass, and is the Bose-Einstein distribution function. The transverse gluon projection 
operator is Vj- = —6^^ + q^q-' /q^, with all other components zero, and T^" are the colour 
SU{?>) matrices. We use the standard high energy conventions that h = c = = 1- 

The integrals over q^ and cos 6 = k ■ q can be evaluated by using the two 5-functions, 
while the remaining angular integrals are all trivial: 

pq 9^ r ^ ( ^ /■^"'-'"' , ^ 

Krp = -— / aujnB{uj) / dqq 



67r4 J2m "'Wo [uj^ - g2 _ RenT(g)] + Au^-f^ 



X 



I dE Tr [(^ + m)7^(^- Q + m)^,] P^- (5) 



The trace over colour indices gives a factor 4/3, the quark energy k'^ = E = (k'^ + m'^Y^'^, 
and the gluon energy g° = u. The transverse component of the trace is 

Tr [(^ + m)7^(^- ^ + m)^,] = 4 piT^iT. - K^Q, - Q^K, + K ■ Qg,,] 

= 4[cu2-g2_2p(l-cos2^)] , (6) 

where 

cos^ = {2ujE-u'^ + q^)/{2kq). (7) 
We then integrate the trace over the allowed range of E, where \uj'^ — q^ — 2Euj\ < 2kq: 



q Ag'^ f°° j-yJuj-^-Am? "frp (cu^ — q"^ + 2rn?) _ 
^T = TT^ duJ UJUBiuj) / dq q —— ^ ^ ^X^,2 "o^ 9 • (8) 
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97r4 72m ' ' Jo [cj2 _ ^2 _ ReHTlQ)] + 4cj2^2 

The integral in eq. (||) is cut off by when u/T ^ 1. The crucial point of our 
analysis is that the leading contribution comes when u and q are hard, of order T. The 
soft Q contribution is perfectly regular and subleading by 0{as)- Because the dominant 
gluon momenta are hard, we make the subsequent approximations: i) jt is independent 
of the momentum and is given by 



IT = ^ In + 1.09681 ... when rrimag > 7t, 

= ^^\n— + OiasTy, (9) 
2 as 

ii) as we are only interested in the leading order, we take Rellr = 3co'q/2, as changing the 
5-function 6{Q'^ — Ren-r) by a Lorentzian with width changes Rellr = 3u!q/2 only by 
0{g'^T'^) corrections. The gluon damping rate depends on the magnetic mass which is, 
for dimensional reasons, a number times g^T |jl3|. Fortunately, this uncertainty enters 
only at the logarithmic level. 
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When m = 0, we evaluate the integral over q analytically, obtaining a simple result: 



T = / duouj nB[uJ) < m — -. ^ — - H arctan h arctan — 



OcUq + 167|.co'2 277-0; y 47^0; 7t 

(10) 

where terms of order higher than have been dropped. Notice that the limit 7t — > 
reproduces the result using a gluon propagator without a finite damping rate. 
The final result is 



.3 = iiM„i A„ ±\ 



It is in fact the dependence of the rate that makes it so large. Should we use a 
bare gluon propagator, S{Q^), the rate would identically vanish. With just the hard 
thermal mass, (5(Q^ — 3^75/2), the rate is of order g^T'^ with no logarithmic dependence. 
Taking into account the anomalously large damping rate 77- shifts the gluon on-shellness 
~ by a logarithmic correction. The additional logarithm has a kinematic origin. 
The processes that lead to our result are of the kind shown in Fig. 3. The enhancement 
comes both from the soft space-like gluon exchange (a Landau-damping term) and from 
the pole of the gluon propagator that is almost on shell, hence the two logarithms. This 
diagram is not suppressed kinematically. Notice also that using a resummed propagator 
for a hard line is not a new feature. For the calculation of the damping rate itself, one 
needs to use the same propagator |p!2| , p!^ . 



It is easy to see that the longitudinal decay is subleading. The only difference (aside 
from the self-energies) between the longitudinal and transverse modes is a factor of two, 
as there are two transverse modes and only one longitudinal mode. At hard momentum, 
the longitudinal thermal mass is exponentially suppressed, ReHj;, = 0, so we have in the 
case of zero bare quark mass 

However, the longitudinal gluon damping rate is of order g'^T with no log. The quark 
production due to longitudinal gluon decay is therefore of the same order as the neglected 
terms in eq. (|TT|) . 

The discussion for massive quarks is the same through eq. (H). When m S> cjq, 
the rate behaves perturbatively as g'^'jr oc g'^lnl/g, and is therefore larger than other 
contributions. We show in Fig. 4 the production rate as a function of m/T for GP 
{Nf = 0); the curves are very slightly altered for Nf = 2. 

We stress again that our result is the leading order perturbatively (when g ^ 0), 
containing all relevant contributions to order aKlnl/a^)^. The fact that the quark 
gluon plasma is not really a system of weak coupling is a different problem, and should 
in principle be solved by adding higher order contributions to our result. One should 
work in this perturbative picture to make a reliable prediction. 

In any case, there seems to be a huge uncertainty in the quark production rate when 
considering large values of g {g > 1). As the rate is linearly dependent on the damping 



rate, which itself depends on the value of the magnetic mass, which is unknown, it is 
difficult at present to evaluate the effect of higher order terms that are important for 
large g. We believe that the rate can easily vary by an order of magnitude. One should 
bear this in mind when making predictions for ultrarelativistic heavy ion collisions. The 
disadvantage of this uncertainty can be turned into a gain, if one is able to measure the 
thermal charm or bottom quark production. This would give a direct measurement of 
the gluon damping rate, and hence of the magnetic mass. 

Having computed the quark production rate, we also discuss the quark equilibration 
time given by Tq = N^q/Rj^, where Ngq is the quark density at equilibrium. In Fig. 5 
we plot Tq as a function of m/T. In calculating Ngg, we take into account the complete 
mass (bare plus thermal corrections) at hard momentum, 

m| = H — when m -^T, 

2027^2 

m| = when m:^T. (13) 
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We cannot rigorously calculate the thermal mass between these limits, so we use the 
ansatz 



2 2 {2 + e-'^/^)g^T^ 

ml^m^ + . (14) 
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For vanishing quark mass, we obtain 

^ 16aU\l\/asy^~'- ^^^^ 

For comparison, we also give the gluon equilibration time. By detailed balance, the 
thermal gluon formation rate must be equal to the decay rate (twice the damping rate), 
so that 

T-e = ^ -, — (16) 

^ 3aslnl/as ^ ' 

We see that gluons equilibrate faster than even the lightest quarks. This stems from the 
fact that quarks are always produced by pairs, which is a relatively slow process, while 
gluons can be copiously emitted by radiative processes. 

Again, our results strongly depend on the value of the coupling constant. From the 
numbers that we obtain, it seems clear that, in ultrarelativistic heavy ion collisions, only 
the light quarks (m, d, and s) may have time to come to chemical equilibrium, while the 
heavy quarks will not. Even taking g( = 3, rq > 10 for massless quarks, so chemical 
equilibration will proceed very slowly. 

In conclusion, we have calculated the quark production rate from pure glue and quark- 
gluon plasma in equilibrium. The leading contribution comes from the decay of trans- 
verse gluons, which are mainly characterized by an anomalously large damping rate 7^. 
The analogy of this process with the plasmon decay into is very close but mislead- 
ing. In the latter case, the photon damping rate is smaller, 7^ = C(e^Tln(l/e)), so 
the plasmon is really described by a quasi-particle with an effective thermal mass. In 
a QCD plasma, this is clearly not the case, and one must be very careful in computing 



the thermal rates by incorporating the damping rate inside the gluon propagator. Our 
result is another example of a physical quantity that is sensitive to the non-perturbative 
magnetic mass scale [|12| . 



Although our calculation is perturbatively correct, there are large uncertainties when 
making predictions for ultra-relativistic heavy ion collisions for the following reasons: i) 
the rate is linearly dependent on jt, which depends on the unknown magnetic mass, 
and ii) the quark gluon plasma, which could be formed in such collisions is not really a 
system of weak coupling, so higher order terms should be calculated. 

This material is based upon work supported by the North Atlantic Treaty Organization 
under a grant awarded in 1991. 
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Figure captions 



1. The lowest-order "cut" diagram for thermal quark production. 

2. The naive lowest order Feynman diagrams for quark- ant iquark production. 

3. A Feynman diagram contributing to leading order for thermal quark production. 

4. The quark production rate, i?^, in a GP {Nf = 0) as a function of m/T. 

5. The quark chemical equilibration time, rq, in a GP as a function of m/T. 
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